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Abstract 

We calculate perturbatively the pressure of a dilute gas of anyons through second order in the 
anyon coupling constant, as described by Chern-Simons field theory. Near Bose statistics, 
the divergences in the perturbative expansion are exactly cancelled by a two-body 5-function 
potential which is not required near Fermi statistics. To the order considered, we find no 
need for a non-hermitian Hamiltonian. 
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This letter is about a two-dimensional non relativistic system of N particles in the ther- 
modynamic limit, interacting with a Chern-Simons gauge field, and interacting among them- 
selves via a two-body 5-function potential. This system has been shown to obey fractional 
statistics [1]. From a quantum mechanical point of view, a system of anyons can be studied 
by introducing a flux tube at the position of each particle, whose magnitude 9 will deter- 
mine the phase associated with relative particle motion. The value 6* = corresponds to 
the anyons being bosons, whereas 9 — tt corresponds to fermions. In the dilute gas regime, 
when the mean interparticle spacing is much larger than the thermal wavelenght, 

\t — (27r/mT)-^/^, the first correction to the pressure comes from the reduced two-body 
hamiltonian. This can be exactly solved in terms of the Aharonov-Bohm phase shifts [2]. 
Since the full spectrum of the N-anyon system is not known, a computation of the pressure 
based on N-cluster expansion [3] is not possible. 

Nevertheless, full knowledge of the spectrum is not neccesary when the anyon model is 
considered from the point of view of quantum field theory. One needs the thermo dynamical 
potential fl{fi,T, V) = —PV as a power series in the fugacity z = exp(/i/T), which can be 
computed via perturbation theory. This procedure is reliable when the field-theoretical model 
accounts for the fractional statistics of anyons. Thus, the question one would like to ask is 
this: which field-theoretical model can be used to reliably compute higher virial coefficients 
(perturbatively in 9)1 

This question has been considered in a series of papers [4-6], where a two-body non- 
hermitian interaction is introduced in order to handle the divergences due to the statistical 
gauge field. As a consequence, the perturbativc treatment must be carefully done, using 
either a harmonic regulator, or considering the system in a box. The source of the non- 
hermitian interaction arises as follows: In the regular gauge where the statistics is standard, 
the anyonic gas is described by a Hamiltonian where three-body interactions of order 9^ are 
present. These lead to divergences in the perturbative expansion. However, the three-boby 
interactions can be removed from the Hamiltonian by a complex gauge transformation [5], 
leaving a new two-body Hamiltonian with a non-hermitian term of order \9\. Then the 
gauge transformation might be understood as performing a self-adjoint extension [9] on the 
original Hamiltonian, representing an additional interaction at order 9, namely a two-body 
5-function. 
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Another viewpoint is afforded if one consider the anyonic gas as described by a Chern- 
Simons gauge field coupled to bosonic or fermionic matter. In this case, the singular three- 
body interactions are kept and they correspond to the last vertex in Table 1. Even so a 
perturbative expansion in 6 will be reliable if divergences cancel out due a two-body 5- 
function. It is not a priori clear whether both approaches must lead to the same results. 

The purpose of this letter is to calculate the contributions to the pressure to order in 
Chern-Simons field theory. This involves a set of three-loop diagrams, whose computation 
to second order in the fugacity is our main new result. Contributions of order e^^^'^ are 
neglected and will be reported in future work. An obvious prerequisite for the validity of 
such computation is the reproduction of the second virial coefficient to order 6"^. We show 
that Chern-Simons field theory with the inclusion of a repulsive 5-function potential among 
particles is a model which exactly reproduces the correct pressure when the fidutial statistics 
is bosonic. The strcnght of this interaction must be adjusted to cancel UV divergences. 
No additional interaction is required when the fidutial statistics is fermionic. Therefore, 
regardless of whether one works with the non-hermitian Hamiltonian of refs. [4-6] or with 
Chern-Simons field theory, the result for the second virial coefficient is the same. 

We shall consider a spinless non relativistic self-coupled field interacting with a Chern- 
Simons gauge field, with Lagrangian density given by 

>C(i,x) = -;^9taxa + -aoS - ^ (V • a)^ ^iI^^iDqiI^ - ^ jDV'P +// V'V - t 
2k k 2p 2m 4 

Do = dt + iao, D = V-^a (1) 

Here S = V x a is the magnetic field, fi is the chemical potential and p is a gauge fixing 
parameter. The Coulomb gauge to be used refers to the choice p — 0. Finally, a — 0{k,) is 
the strength of the two-body contact interaction of the form 

a 

y(xi-X2) = -5(xi-X2). (2) 

When the starting statistics is fermionic we will set a — 0. 

In the imaginary-time formalism of thermal field theory [7], the functional integral ex- 
pression for the grand partition function involves an integration over imaginary time from 
to (3^ T-i , 



Z{p,T,V)= I VaoVaVij^Vip exp ( f dr I d\ L{t = -iT,x)\ 

J (anti)periodic ^Jo J ' 



(3) 



Here (anti)periodic means that the integration over fields is constrained so that P) — 
±'0(x, 0) where the (lower) upper sign refers to (fermions) bosons. 

Now we can proceed to expand in a power series in n, by using a set of diagrammatic 
rules which follows from (1). These are listed in Table 1. The diagrams describing the 
perturbative series for \nZ have the form of connected closed loops. It should be noted that 
there will be a factor of left over for each graph, corresponding to the extensivity of In Z. 
This factor cancel out in expressions for the pressure. 

The non-zero graphs contributing to the pressure up to order are shown in fig.l. 
Note that all graphs of order k vanish due to the index summatios associated with vertices. 

The same argument applies to graphs of order an, although at higher orders this is not 
necessarily so. In order to keep the correct order of the operators according to their r-vahics, 
one must insert a factor e*"^"^ whenever a particle line either closes on itself or is joined by 
the same instantaneous interaction line. We take r/ — O"*" at the end of the calculation. As 
a consequence, terms involving three bubbles in P^^ and in pi^-* are of order e^i^f^. This 
follows simply from power counting with the formulae 



-C — n^TTTTTT--. — rr^ = -C"c(q). W 



-T 

P'^iun + fi- qV(2m) exp [/3 (qV(2m) - li)] - ( 

-E- 



^Z^i,„ + /q2/(2„)=-^C"<(q). (5) 

where ^ = 1 for bosons and ^ = — 1 for fermions. Therefore, graphs (c) and (g) can be 
neglected, since we are computing to order e^^^. 

We are now in position to derive the perturbative corrections to the pressure 

P = ^lnZ. (6) 



Graph (a) gives 

(1) [ . . . . _ 3/3M^ 

2 7 (27r)2 (27r)2''^^P^''^^'^^~ 2X^^ ^ 



and there are not first order correction with fermionic statistics. 
The contribution from graph (b) may be written 

.,2 r ^2 



p(2) a /" q r-rrO / m2 
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where 

(fp n^(p + q) -n<;(p) 



(27r)2 iz/„-ep+q + £p 



(9) 



represents the lowest-order density-density correlator (see fig.2). At high temperature and 
low density, distribution functions reduce to n^(p) = e^^e"^"^ Z^"*. Then, keeping terms to 
second order in the fugacity, eq.(8) gives 



where = \/ q* is a cutoff and $(a;) is the plasma dispersion function 



p"*^ / 2 2 

$(x)=2/ dt ^2e-^ / die^ (11) 

Jo vx^ — Jo 

^{x) has the following limiting behavior 



(12) 



$(a;) = 2x + 0(x"^) x < 1, 
$(x) = x"^ 0(x"^) a;>l. 

We have found a UV logarithmic divergence. This corresponds to the divergence in the 
second Born approximation for the scattering amplitude with a 5-function potential [9] . 

Now we consider the ring contribution from graph (d) of fig. 1, 

e2 ^ /• d\ nO(i/n, q) E^iun, q) 



.d_ /• ^'q n>n,q)sVn,q) 

" 2(3^ J (27r)2 q2 ' ^^''^ 



where 'EP{i'n, q) is the two-dimensional transverse component of the lowest-order current- 
current correlation given by (see fig.2) 

n%-(^n,q)=-^^E /7^[(P+f)z(P+f),-^"(^l,p)^°(^l + ^n,q+p) 

= -n%o (-n, q) 4 ^ + -° (-n, q) (5,,- - ^) . (14) 



In the classical limit, when n(^{p) — e^^e /^P^/2"i^ 
and eq. (13) gives 



(15) 



pa 



K e 



r 


"$(a;) 2" 


/ dx 




Jo 


x-^ a; 



(16) 



Again a UV logarithmic divergence appears coming from the ring loop with a seagull vertex. 
We still have to compute the exchange graphs (e) and (f) of fig. 1. These give 

(Pq (Pp d'^k 



P. = -c 



K 



E 



4m2/33 ^ J (27r)2 (27r)2 (27r)2 



6?°(a;i,p) 6?°(u;i + i/,q + p) 



X ^0(a;2 - k - q) ^0(o.2, k) ^ ^^^^.^ ^/ , 

(p - k + q)^ q^ 



(17) 



K 



E / 7^ 7^ 7^ P) ^°(^2, k) 



^(^\trJ ^27r)2 (27r)2(27r)2 



X ^?°(a;i + t/,p + q) 



(P - k + q) ■ q 



(p - k + q)2 q2 ■ 

A detailed evaluation of these terms to order e^^^ yields the divergent contribution 

2 2/3^ rx* 



(18) 



' 47rm Ai Jo 



(19) 



where the divergence comes again from the graph (f) with a seagull vertex. 

In the fermionic case, the total contribution from the gauge field is finite by itself, giving 



p(2) — pd , pe , pf 



At2e2/3M 

Snm Ay Jq 
STTm Ay 

^2e2/^M 

Anm Ay 



CO 


"$(a; 


) 2 


dx 








a;2 


a; 






"$(a;)" 


1 dx — 


JQ 







+ 2$(a;) 



(20) 



In the bosonic case, the divergences in the total contribution to the pressure to order 
cancel out only if a = -^2K/m, and we find 



(2) 



dx 



' ' - -2$(x)+4$(x) 



8nm Ay Jq \_ x^ x 
iirm Ay ' 

which is the same result as in fermionic case. 

Finally, using the pressure for two-dimensional ideal quantum gases 

pW(/.,r) = CrAy2Li2(Ce'^/^), 

where IA2 denotes the dilogarithm function, we obtain the cluster expansion 
P(/i,T) 1 V 



(21) 



(22) 



^Y.^iz'= T-HP(°)(/., T) + P^'^ + + Oi 

1=1 



K 



A^ An 



:i+c)«^ « 



47r 



87r2 



(23) 



When the fidutial statistics is bosonic this formula reproduces the correct second virial coef- 
ficient by taking the upper sign corresponding to a repulsive contact interaction of strength 
a — 2k /m and putting k — 26. The other possible choice, a — — 2«;/m, corresponds to the 
nonrelativistic model admitting classical self-dual solutions [10]. When the fidutial statistics 
is fermionic, the suitable choice is k — 2{9 — tt), where 9 is the statistical parameter in both 
cases. 

To conclude, we have shown that standard perturbative Chern-Simons field theory cou- 
pled to nonrelativistic matter exactly reproduces the second virial coefficient. Through 
second order in 0, all divergences cancel, and there is no need for renormalization or re- 
summation. A non-hermitian Hamiltonian is not required to account for the statistical 
interaction between anyons. 
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Ta.l)le 1. Dare 'prvpacjators and vertices 



^?°(a;n,p) = -. ■ ^ 0//0 ^ 

2m: . , 
cun — —5— m boson propagator 

(2n+l)7r . 
cUfi — m rermion propagator 



^00 (q) = =0 in the Coulomb gauge, p = 



r" = -a 



ro = -i 



Pj + kj 



_ J_ 
m 
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Figure Captions 

Fig. 1. Non zero diagrams contributing to the pressure to the second order in k. The sign 
± refers to Bose or Fermi propagators. Graphs (c) and (g) contribute to the third virial 
coefficient. Combinatoric factors are shown in the diagram. 

Fig. 2. The self-energy of the statistical gauge field at the one- loop level in terms of the 
functional derivative of the (zero) pressure to lowest order. IPI means that only the one- 
particle irreducible diagrams contribute to liP. 
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